Pauli susceptibility of A3C60 (A=K, Rb) 
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The Pauli paramagnetic susceptibility of A3C60 (A= K, Rb) compounds is calculated. A lattice 
quantum Monte Carlo method is applied to a multi-band Hubbard model, including the on-site 
Coulomb interaction U . It is found that the many-body enhancement of the susceptibility is of the 
order of a factor of three. This reconciles estimates of the density of states from the susceptibility 
with other estimates. The enhancement is an example of a substantial many-body effect in the 
doped fullerenes. 
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The Pauli paramagnetic susceptibility x is interesting 
for several reasons. Firstly, x is enhanced by many-body 
effects relative to its value xo for noninteracting elec- 
trons, and x/xo is one measure of the strength of the 
many-body effects in the system. In the alkali-doped C@o 
compounds, A3C60 (A= K, Rb) the Coulomb interaction 
U between two electrons on the same molecule is largetl 
compared with the width W of the partlv. filled f i„ band, 
with typical estimates U/W ~ 1.5 — 2.5.H In view of this 
large ratio, one expects very strong many-body effects 
for these systems. Up to now, however, there seems to 
be no unambiguous signature of such strong effects. 

Secondly, xo is related to the density of states N(0) at 
the Fermi energy, and values of N(0) can be extracted 
from x if the enhancement x/xo is known. A^(0) is 
important for the superconductivity and the electron- 
phonon interaction A, since theoretical calculations give 
X/N(0), while some experimental (e.g., neutronEl and Ra- 
man scatteringa) estimates give AAi(O) and others (e.g., 
photoemission from C^~ moleculesB) give X/N(0). The 
estimate of N(0) is therefore crucial for obtaining val- 
ues of A and for our understanding of the supercon- 
ductivity. Typically, for K3C60 band structure calcula- 
tions give N(0) ^6 — 9 states/ (eV-spin)lrO Estimates 
based on the specific heat and the NMR relaxation rate 
give N(0) ~ 5 -j£ states/ (eV-spin)B and N{0) ~ 7.2 
states/ (eV-spin),E3 respectively. On the other hand, 
much larger values N(Q) ~ 10 — 16 statesZ(£i£|Spin) for 
K3C60 are deduced from the susceptibility,E-3tZl if many- 
body effects are neglected. A substantial many;-body en- 
hancement (factor 2-3) for the susceptibility,£3 could es- 
sentially reconcile these rather different estimates. On 
the other hand, density functional calculations in the lo- 
cal density approximation (LDA) find that the enhance- 
ment is only about a factor of 1.3-1.4.E3 We note that 
N(0) here refers to densities of states defined in some- 
what different ways depending on the experiment, as dis- 
cussed below. 

We tewe used a lattice quantum Monte Carlo 
methodtSH for calculating the Pauli susceptibility for a 



multi-band Hubbard model of the system. For a small 
system with four Cgo molecules, we demonstrate that this 
method gives an accurate enhancement of the suscepti- 
bility. For realistic values of the parameters, the sus- 
ceptibility is enhanced by about a factor of three, which 
essentially reconciles estimates of N(0) based on the sus- 
ceptibility with other estimates. 

In the presence of a small external magnetic field H, 
the energy of the system can be written as 

E{M) = E (M) - MH « £00 + J«M 2 - MH, (1) 

where M = — //s(iVf — N±) is the magnetic moment of 
the system, with [Ib being the Bohr magneton and N a 
the number of electrons with spin a. Eq(M) is the en- 
ergy of the system with a moment M in the absence of 
an external field. Minimizing the energy with respect to 
M, we obtain the susceptibility 
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X ~ H a 



(2) 



In the following we therefore calculate Eq(M) for the 
interacting and noninteracting (U — 0) systems, from 
which we obtain the many-body enhancement x/xo- 

We use a multi-band Hubbard model of the A 3 Ceo 
compounds 
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where the first term describes the three-fold degener- 
ate t\ u states on the sites i and with orbital (m) and 
spin (a) indices. The second term describes the hop- 
ping between the sites, and the third term describes 
the Coulomb on-site interaction. Multiplet effects and 
the electron-phonon interaction have been neglected. In 
A3C60, the Ceo molecules are preferentially in one pf,two 
possible orientations in an essentially random way.t3 We 
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take this into account by having a large cell, where each 
molecule takes one of the two preferred orientations in a 
random way, and the hopping matrix elements between 
two molecules itate into account the orientations of these 
two molecules 

To calculate the energy of the model in Eq. (|^), we 
use a T = projection lattice Monte Carlo method, 
introduced by ten Haaf et alHH In this method a trial 
function is constructed from a Slater determinant, using 
a Gutzwiller Ansatz.Ea An approximate ground-state is 
then projected out in a diffusion Monte Carlo (DMC) 
approach, using a "fixed node" approximation. This 
method has been used to stud* the condition for a Mott- 
Hubbard transition in A 3 C6o-n 

To test the accuracy of the DMC approach, we 
have first applied the method to a cluster of four Ceo 
molecules. This cluster is so small that we can also ob- 
tain the exact solution for the model in Eq. (|^), using 
exact diagonalization. We then calculate the coefficient 
a in Eq. (Q) by considering the energy for iV-p — N± — 
and 2. The DMC and exact results are compared in 
Table Q. We can see that the DMC method is quite ac- 
curate in this case, and if a similar accuracy is obtained 
for larger systems, it is quite sufficient. 

TABLE I. The enhancement of the susceptibility for a 
model with four Ceo molecules according to diffusion Monte 
Carlo (DMC) and exact calculations as a function of the 
Coulomb energy U. The band width is W = 0.58 eV. 
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similar results have been obtained for a Hubbard model 
without orbital degeneracy and in the limit of infinite 
dimensions .c3 
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FIG. 1. The energy as a function of iVf — ~ M rela- 
tive to the energy for M=0. Second order curves have been 
fitted to the energies. Results for U = 0, 0.5, 1.0, 1.5 and 
2.0 eV are shown. The calculations are for 32 molecules and 
the energy is in eV. 



Within the Hartree-Fock approximation, the suscepti- 
bility behaves as 



X 3 



The enhancement of the susceptibility is sensitive to 
the density of states (DOS) close to the Fermi energy Ep . 
For small and intermediate size clusters of Cqq molecules, 
the DOS depends on the orientations of the molecules, 
while for large clusters the DOS rapidly converges. Since 
we can only treat intermediate size clusters (~ (32 — 64) 
molecules) in DMC, we have therefore chosen orienta- 
tions which in a one-particle approximation give similar 
DOS close to Ep as for very large clusters. 

In Fig. [I] we show results for the total energy as a 
function of M ~ Nf — JVj_ for different values of U. The 
results can be rather well fitted by parabolas, although 
the precise parameters of the parabolas have a certain 
dependence on the range of M considered. From these 
slopes we can immediately deduce values of the enhance- 
ment x/Xo- I n Fig- H the inverse of the enhancement 
Xo/x is shown. It is immediately clear that the enhance- 
ment grows with U and that x would diverge for U a 
bit larger than 2 eV, if no other transition (e.g., antifer- 
romagnetic) happened before. Estimates of U are typi- 
cally in the range 1-1.5 eV, giving an enhancement of the 
susceptibility by about a factor of three. Qualitatively 



where the factor three comes from the three-fold degen- 
eracy and N(0) is the density of states per spin, which 
here is about 5.5 states/ (eV-spin). Here we have assumed 
that the three orbitals are equivalent. The DMC results 
show a similar behavior, but with a prefactor in front 
of U , which is between a factor of four and five times 
smaller. This large change in the prefactor illustrates 
the importance of correlation in these systems. 

To deduce N(0) we need to know the Pauli (param- 
agnetic) susceptibility. Measurements using a SQUID 
may also contain a diamagnetic contribution, while EPR 
measurements do not. In Table || we show various exper- 
imental results converted to N(0). From the SQUID re- 
sults, diamagnetic contributions estimated by the respec- 
tive authors have been subtracted, but the many-body 
enhancement has not been considered. We can see that 
the results range between 10 and 16 states/ (eV-spin). 
If we consider a many-body enhancement of a factor of 
three, as deduced above, these results would be reduced 
to about N susc (0) ~ 4 - 5 states/ (eV-spin). 
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FIG. 2. The inverse enhancement Xo/x °f the suscepti- 
bility for a cluster with 32 molecules as a function of the 
Coulomb interaction U (in eV). The band width is 0.66 eV. 

Up to now we have not fully considered differences in 
the definition of the density of states N(0) deduced from 
different sources. Calculations of the electron-phonon in- 
teraction A are usually based on band structure calcu- 
lations, giving a density of states N (0) for noninteract- 
ing electrons. The electron-electron interaction modifies 
iVo(O) to its interacting value 7V(0) (density of quasi- 
particle states), which should be used in the calculation 
of A.E3 This density of states enters in the specific heatcll 



A/-(0)(1 + A), 



(5) 



where the last factor is due to the electron-phonon inter- 
action, assuming that the hand width is large compared 
with the phonon frequencyHj In the-same approximation, 
the susceptibility can be written aaHJ 



where F§ is a Landau parameter and 1/(1 + Fft) can 
be considered as a Stoner enhancement. Our calcula- 
tion x/chio includes both the factors A/"(0)/Ao(0) and 
(1 + Fq). Within the present Monto Carlo technique 
it is not possible to calculate the specific heat, and we 
can therefore not separate the two contributions. The 
value of N SUSC (Q) deduced from the susceptibility, after 
dividing out many-body effects, should therefore primar- 
ily be compared with the results obtained from band 
structure calculations. The band structure calculations 
have been performed for orientationally ordered systems 
and give Nq 1 " 4 ~ 6 — 9 states/ (eV-spin). Sine the real 
systems have orientational disorder, we estimate the cor- 
responding density of states by solving the Hamiltonian 
(||) for U = and orientational disorder. Comparision 
with calculations for ordered systems suggests a reduc- 
tion of N(Q) by about 15 - 20% to Nf sord (0) -5-7 
states/ (eV-spin). This is in rather good agreement with 
N susc (0) - 4 - 5 states/ (eV-spin). 

It is interesting that the estimate oOZ (0) ~ 5 — 6 states 
/(eV-spin) based on the specific heatJla is comparable to 



the noninteracting result N$ isord (0) - 5 - 7 states/ (eV- 
spin). In contrast to what has been found for the-.non- 
degenerate Hubbard model in infinite dimensions,!^ this 
surprising observation suggests that the enhancement of 
the density of states is small in the A3C60 compounds or 
that there may even be a reduction. If this is indeed the 
case, then the enhancement of the susceptibility would be 
a Stoner enhancement, and the Nq usc (Q) obtained from 
the susceptibility, after dividing out the many-body en- 
hancement, could be compared with other experimen- 
tal estimates of iV(0). We then find that these exper- 
imental estimates are essentially brought in line with 
each other and with the band structure estimates, giving 
N(0) - 5 - 7 states/ (eV-spin) for K 3 C 60 . 

In our model and in the calculations of the sus- 
ceptibility, we have neglected multiplct effects and the 
electron-phonon interaction, which are now discussed. 
The electron-phonon interaction leads to an increase of 
the density of states at the Fermi energy, due to the re- 
duced dispersion of states within roughly a phonon en- 
ergy of the Fermi energy. This does not influence the 
susceptibility, if the phonon-energies are small compared 
with the electronic energiesO Although this assumption 
may not be entirely satisfied for A3C60 and interesting ef- 
fects may result from the finite band width, c3 we here ne- 
glect the effects of the electron-phonon interaction on the 
density of states. Instead we focus on how the electron- 
phonon interaction influences the moment formation on 
the Ceo molecules. 

TABLE II. The density of states N(Q) (per eV and spin) 
for K3C60 as deduced from susceptibility measurements. The 
results have not been corrected for the Stoner enhancement, 
which would lead to reduced estimates of N(0). 
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If the multiplets are neglected but the electron-phonon 
interaction (Jahn- Teller effect) is considered, the lowest 
spin 1/2 state is favored over the spin 3/2 state, accord- 
ing to calculations £ar the lowest state of each multiplicity 
for a free moleculeE3 The energy lowering of the spin 1/2 
relative to the spin 3/2 state may be as large as 0.3 eV.EHl 
This suppresses the formation of moments and probably 
tends to reduce the susceptibility. On the other hand, the 
multiplet effects should favor the formation of moments 
on the Ceo molecules, by giving preference to states with 
the spin 3/2. 

The multiplet effects lead to five spin 1/2 states with 
the energy 3K and three spin 1/2 states with the en- 
ergy 5K relative to the spin 3/2 state. Here K is the 
exchange integral between two t\ u orbitals, and we have 
assumed that the Coulomb integral (U xx ) between equal 
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orbitals is 2K larger than the one (U xy ) between different 
orbitals. To estimate K, we have used a simple model, 
where the Coulomb integral between two 2p-charge dis- 
tributions on two carbon atoms goes as A'-R, where R 
is the separation between the two atoms.Eil The on-site 
interaction was assumed to be 15 eV. Without screen- 
ing we find that K = 0.12 eV. This probably overesti- 
mates the multiplet effects. An alternative estimate is 
obtained by using RPA screening of the Coulomb inter- 
action. We then find K — 0.030 eV. A similar result 
(0.024 eV) was alsOripund by Joubert, using a density 
functional approach H3 

From the numbers above it follows that there should 
be a partial cancellation between electron-phonon and 
multiplet effects. Depending on which numbers are used, 
either effect could be argued to be larger. If the electron- 
phonon effects win, this may lead to a somewhat smaller 
enhancement of the susceptibility than was found above 
(Fig. 1). 

We have calculated the Pauli susceptibility of the 
doped fullerenes A 3 C 6 o (A= K, Rb). The enhancement 
is of the order of a factor of three, which allows us to 
reconcile the estimates of the density of states from the 
susceptibility with other estimates. This suggests that 
for K 3 C 60 N(0) ~ 5 - 7 states/ (eV-spin). This value 
is only slightly smaller than a value (A^(0) = 7.2) used 
recently to provide support for an electron-phoxion mech- 
anism driving the superconductivity in K3C60P but sub- 
stantially smaller than some values used in early the- 
oretical discussions. The susceptibility enhancement is 
appreciably larger than the one (factor 1.3-1.4) found in 
the LDA, and it is one of the first explicit examples of 
important many-body effects, expected to be found in 
these systems. Comparison with Hartree-Fock calcula- 
tions shows, however, that the enhancement is about four 
to five times smaller than the HF result, illustrating the 
importance of correlation effects. 
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